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EXECUTIVE  SUMMARY 


Water  droplet  ingestion  into  turbine  engines  resulting  from  heavy  rain  and 
wheel  spray  generated  on  a  wet  runway  is  of  importance.  Adverse  effects  of 
large  quantities  of  water  ingestion  can  include  the  compressor  stall  and 
combustor  flame-out.  Engine  certification  requirements  as  set  forth  in  FAA 
regulations  call  for  continued  engine  operation  at  takeoff  and  flight  idle 
conditions  while  ingesting  water  at  4  percent  by  weight  of  airflow,  generated 
by  a  spray  to  simulate  rain.  There  is  also  a  certification  requirement  on  the 
entire  aircraft  system  which  dictates  that  the  system  must  be  designed  to 
prevent  hazardous  quantities  of  water  from  being  ingested  into  the  engine 
during  takeoff,  landing  and  taxiing  operations  on  wet  runways.  The  present 
work  was  undertaken  to  develop  measurement  techniques  of  two-phase  droplet 
laden  airstreams  during  engine  water  ingestion.  The  ultimate  objective  is  to 
correlate  the  non- intrusive  measurements  of  the  water  ingestion  rate  and 
droplet  size  and  spatial  distribution  at  the  engine  inlet  with  engine 
performance  parameters.  Such  techniques  and  data  will  assist  the  FAA  in 
evaluating  current  water  ingestion  certification  tests. 

A  r.on- intrus ive  optical  technique  was  developed  for  the  determination  of  liquid 
mass  flux  in  a  droplet  laden  airstream.  The  technique  is  also  capable  of 
providing  information  on  the  droplet  size  and  spatial  distribution  at  the 
nacelle  inlet  plane. 

Independent  measurements  of  the  liquid  water  content  (LWC)  of  the  droplet  laden 
airstream  and  of  the  droplet  velocities  were  made  at  the  inlet  plane  of  a 
simulated  nacelle  in  a  wind  tunnel  for  the  liquid  mass  flux  determination.  The 
liquid  water  content  was  determined  by  illuminating  and  photographing  the 
droplets  contained  within  a  thin  slice  of  the  flow-field  by  means  of  a  sheet  of 
light  from  a  pulsed  laser.  Fluorescent  dye  introduced  in  the  water  enhanced 
the  droplet  image  definition.  The  droplet  velocities  were  determined  from 
double  exposed  photographs  of  the  moving  droplet  field.  The  technique  was 
initially  applied  to  a  steady  spray  generated  in  a  wind  tunnel.  It  was  found 
that  although  the  spray  was  initially  steady,  the  aerodynamic  breakup  process 
was  inherently  unsteady.  This  resulted  in  a  wide  variation  of  the  instantane¬ 
ous  liquid  water  content  of  the  droplet  laden  airstream.  The  standard 
deviation  of  ten  separate  LWC  measurements  was  31  percent  of  the  average. 
However,  the  liquid  mass  flux  calculated  from  the  average  LWC  and  droplet 
velocities  came  within  10  percent  of  the  Known  water  ingestion  rate. 


INTRODUCTION 


The  effects  of  water  droplet  ingestion  into  turbine  engines  resulting  from 
heavy  rain  and  wheel  spray  generated  on  a  wet  runway  is  of  importance.  The 
effects  of  water  ingestion  on  engine  performance  have  recently  been 
investigated  (reference  1)  and  a  probe  for  stagnation  pressure  measurement  in  a 
droplet  laden  airflow  was  developed  (reference  2)  .  The  adverse  effects  of 
large  quantities  of  water  ingestion  can  include  the  compressor  stall  and 
combustor  flame-out  (reference  1).  Engine  certification  requirement  as  set 
forth  in  FAA  regulations  (reference  3)  call  for  continued  engine  operation  at 
takeoff  and  flight  idle  conditions  while  ingesting  water  at  4  percent  by  weight 
of  airflow,  generated  by  a  spray  to  simulate  rain.  There  is  also  a 
certification  requirement  on  the  entire  aircraft  system  which  dictates  that  the 
system  must  be  designed  to  prevent  hazardous  quantities  of  water  from  being 
ingested  into  the  engine  during  takeoff,  landing,  and  taxiing  operations, 
(reference  4) . 

The  criteria  for  both  water  ingestion  certification  tests  are  somewhat 
arbitrary.  The  4  percent  by  weight  water  ingestion  test  for  engine 
certification  does  not  address  such  issues  as  droplet  size  and  their  spatial 
distribution  over  the  frontal  area  of  the  nacelle.  For  performance  on  a  wet 
runway,  the  criterion  is  even  more  arbitrary  in  that  it  does  not  specify  what 
constitutes  a  hazardous  quantity  of  water  ingestion. 

The  present  work  was  undertaken  to  develop  measurement  techniques  in  two-phase 
droplet  laden  airstreams  to  better  quantify  the  engine  water  ingestion.  The 
ultimate  objective  is  to  correlate  non- intrusive  measurements  of  the  water 
ingestion  rate  and  droplet  size  and  spatial  distribution  at  the  engine  inlet 
with  engine  performance  parameters.  Such  techniques  and  data  will  assist  the 
FAA  in  evaluating  current  water  ingestion  certification  tests. 


REVIEW  OF  APPLICABLE  MEASUREMENT  TECHNIQUES 

A  review  of  available  measurement  techniques  for  sprays  was  undertaken  to 
determine  the  applicability  of  the  existing  techniques  to  quantitative 
determination  of  water  ingestion  into  an  engine.  The  quantities  that  need  to 
be  determined  are: 

a)  Drop  size  distribution  together  with  spatial  distribution  of  drops  at 
the  nacelle  inlet. 

b)  The  mass  flow  rate  of  water  crossing  the  nacelle  inlet  plane  at  any 
instant . 

Extensive  work  on  measurements  of  wheel  sprays  generated  by  aircraft 
under-carriages  was  carried  out  by  Barrett  (reference  5).  A  spray  intensity 
probe  was  developed  to  measure  mean  local  dynamic  pressure  generated  by  moving 
droplets.  However,  all  of  Barrett's  measurements  were  in  the  near  field  of  the 
wheel  generating  the  spray  and  the  technique  used  only  provided  time  averaged 
local  measu  •".ents  In  the  spray. 

Several  mechanical,  electrical,  and  optical  methods  are  available  for  droplet 
size  determination  in  fuel  sprays  as  surveyed  in  review  articles  by  Jones 


(reference  6)  and,  McCreath  and  Beer  (reference  7).  None  of  the  techniques 
surveyed  would  be  capable,  in  their  existing  form,  of  satisfying  the  second 
requirement  above.  It  was  considered  possible  that  the  requirement  could  be 
met  by  some  modification  of  the  existing  techniques. 

The  first  trial  approach  was  an  extension  of  the  charged  wire  probe  technique 
described  by  Gardiner  (reference  8)  for  drop  size  determination  in  water 
sprays.  In  that  technique,  the  electrical  pulse  generated  in  a  circuit 
containing  a  charged  electrode  when  a  drop  impacts  the  electrode  is  measured  by 
a  pulse  height  analyzer.  The  probe  is  initially  calibrated  using  known  size 
drops  impacting  the  electrode.  A  relationship  is  established  between  the  pulse 
height  and  the  drop  size.  In  practice,  the  drop  size  characteristics  are 
derived  from  the  pulse  height  statistics  stored  in  a  pulse  height  analyzer. 
The  technique  as  described  by  Gardiner  (reference  8)  thus  provides  the  local 
drop  size  distribution  in  a  spray.  In  the  present  work,  a  modification  of  the 
technique  was  considered.  Instead  of  a  single  electrode  probe,  a  charged  grid 
was  considered.  A  copper  wire  mesh  with  wire  spacing  of  approximately  1  mm  was 
charged  to  2000  volts  by  a  high  voltage  DC  power  supply.  The  idea  was  to 
install  the  screen  at  the  nacelle  inlet  plane  so  that  all  drops  larger  than  the 
mesh  size  would  be  intercepted  by  the  screen.  The  charge  transfer  between  the 
screen  and  the  impacting  drops  would  set  up  a  current  in  the  circuit  supplying 
the  grid.  The  current  would  be  proportional  to  the  total  surface  area  of  the 
drops  impacting  per  unit  time.  Then  if  the  size  distribution  were  to  be 
determined  by  an  independent  optical  technique,  the  volume  flow  rate  of  water 
would  be  proportional  to  the  product  of  the  Sauter  Mean  Diameter  (SMD)  of  the 
spray  and  the  current  supplied  to  the  grid. 

In  practice,  however,  several  difficulties  were  encountered  with  this 
technique . 

1)  The  charge  distribution  on  the  screen  was  non-uniform,  resulting  in 
different  pulse  characteristics  for  different  impact  locations  on  the 
screen  for  the  same  size  drop. 

2)  Wetting  of  the  screen  altered  the  initial  charge  distribution. 

1)  Water  film  on  the  support  set  up  a  conduction  path  from  the  screen  to 

the  ground,  causing  a  leakage  current. 

Because  of  these  difficulties  and  the  fact  that  an  optical  technique  was  still 
needl'd  in  conjunction  with  the  charged  screen  for  water  flow  determination, 
this  approach  was  abandoned  in  favor  of  a  purely  optical  non- intrus ive 

t  echn i que . 


PRESENT  TECHNIQUE 

A  non- intrusive  optical  technique  was  developed  for  the  determination  of  drop 
sizes,  spatial  distribution  of  drops  at  the  nacelle  inlet  plane  and 
instantaneous  mass  flow  rate  of  liquid  water  entering  the  nacelle.  The  liquid 
water  mass  flow  rate  is  determined  by  independent  measurements  of  the  liquid 
water  content  of  the  droplet  laden  airstream  and  the  droplet  velocity  at  the 
nacelle  inlet  plane. 


For  liquid  water  content  (LWC)  determination,  a  thin  cross-section  of  the  flow 
close  to  the  nacelle  inlet  plane  is  illuminated  by  a  pulsed  laser  light  sheet. 
The  drops  contained  within  this  light  sheet  are  photographed  by  a  camera  placed 
with  its  axis  nearly  normal  to  the  plane  of  the  light  sheet  (figure  1).  The 
duration  of  illumination  of  drops  by  the  laser  pulse  is  extremely  short,  about 
10  nano  second.  Therefore,  the  motion  of  the  drops  is  frozen  in  the 
photographs.  The  thickness  of  the  laser  sheet  is  controlled  by  a  beam  expander 
and  a  slit.  The  slit  was  of  9  mm  width  while  the  initial  beam  width  was  1'*  mm. 
The  slit  thus  allowed  only  the  intense  central  portion  of  the  expanded  beam  to 
pass  through,  cutting  off  the  less  intense  outer  portions. 

The  depth  of  field  of  the  camera  is  set  to  be  larger  than  the  light  sheet 
thickness  by  proper  selection  of  aperture  and  magnification.  This  assures  that 
all  illuminated  drops  (contained  within  the  light  sheet)  appear  in  the 
photographs  with,  a  sharp  edge  definition.  The  measurement  volume  is  defined  by 
the  product  of  the  orojected  frontal  area  of  the  nacelle  and  the  light  sheet 
thickness.  Information  on  drop  sizes  and  spatial  distribution  is  also  obtained 
from  this  photograph. 

•''her.  drops  are  photographed  using  the  optical  set  up  as  shown  in  figure  1,  with 
proper  focusing  and  depth  of  field  selection,  the  droplet  edge  definitions  in 
the  resulting  photographic  images  are  still  far  from  ideal  as  shown  in 
figure  2.  This  is  caused  by  two  effects: 

1)  The  illumination  of  individual  drops  is  not  uniform,  resulting  in 
poor  image  definition  of  larger  drops. 

2)  Scattering  and  diffraction  around  smaller  drops  causes  a  "halo" 
effect  around  the  drop  boundary,  causing  the  image  to  appear  much 
larger  than  the  actual  size. 

An  improved  version  of  this  technique  utilizes  laser  induced  fluorescence  of  a 
small  quantity  (less  than  10  ppm)  of  dye  (Rhodamine  6G  or  Fluorescene) 
introduced  in  the  water.  The  fluorescence  spectrum  of  the  dye  lies  in  a 
wavelength  range  longer  than  that  of  the  incident  light,  therefore,  the 
incident  light  scattered  from  the  drops  can  be  filtered  out  and  only  the 
fluoresced  light  is  photographed.  Shott  colored  glass  filters  are  ideally 
suited  for  this  purpose.  This  technique  results  in  much  better  edge  definition 
of  droplets  in  the  photographic  images  and  a  nearly  uniform  illumination  of 
droplet  interior  (see  figure  3). 

An  automated  digital  image  processing  system  was  developed  to  analyze  the 
photographic  images.  The  image  processing  algorithm  detected  drop  edges, 
defined  droplet  boundaries,  calculated  the  area  of  the  drop  images  and  an 
equivalent  drop  diameter  based  on  this  area.  Droplet  size  distributions  were 
constructed  once  statistics  on  a  sufficient  number  of  drops  were  available. 
The  liquid  water  content  was  determined  from  the  total  volume  of  drops 
contained  within  the  projected  area  of  the  nacelle,  together  with  the  knowledge 
of  the  incident  light  sheet  thickness.  A  description  of  the  image  processing 
system  is  presented  later  in  this  report. 

For  the  determination  of  instantaneous  mass  flow  rate  of  airborne  liquid  water, 
a  measurement  of  drop  velocities  is  needed  in  addition  to  the  LWC  measurement. 
This  may  be  accomplished  non- intrusively  by  laser  double  pulse  photography  of 
the  moving  droplet  field.  Two  images  of  each  drop  appear  in  the  photograph 
(see  figure  A)  and  the  drop  velocity  may  be  determined  by  the  measurement  of 
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the  distance  translated  in  a  known  time  interval.  Solid-state  lasers  such  as 
Ruby  or  Neodymium:  yttrium  aluminum  garnet  (Nd:YAG)  lasers  may  be  operated  in 
the  double  pulse  mode,  wherein  the  Q- switch  is  opened  twice  in  rapid  succession 
to  split  the  energy  of  the  flash  tube  between  two  intense  pulses.  The  time 
interval  between  the  two  pulses  may  be  adjusted  in  the  range  of  50  to  200  ^s . 

The  product  of  the  LWC  evaluated  over  the  nacelle  frontal  area,  the  drop 
velocity  and  the  nacelle  frontal  area  is  the  instantaneous  mass  flow  rate  of 
airborne  liquid  water  into  the  nacelle. 


WATER  INGESTION  SIMULATION  FACILITY 

In  the  first  phase  of  the  program,  an  experimental  facility  was  developed  to 
simulate  engine  water  ingestion  and  to  calibrate  the  non- intrusive  optical 
technique  for  the  quantitative  measurements  of  engine  water  ingestion.  A 
description  of  the  facility  follows. 

The  engine  water  ingestion  simulation  was  set  up  in  an  open  circuit  wind  tunnel 
of  18-inch  x  18-inch  test  section.  A  schematic  diagram  of  this  facility  is 
shown  in  figure  5.  Air  was  supplied  to  the  test  section  by  a  high  capacity 
blower  via  a  settling  chamber,  screens  and  a  nozzle  contraction.  The  blower 
output  was  adjusted  by  a  damper  vane  ring  at  the  inlet.  A  water  spray  nozzle 
was  mounted  on  a  streamlined  sting  support  at  the  entrance  to  the  test  section. 
The  nozzle  consisted  of  seven  2  cm  long  tubes  of  1.6  mm  inside  diameter.  The 
water  was  supplied  to  the  nozzle  from  a  pressurized  tank  via  a  flowmeter  and  a 
ball  valve.  The  individual  tubes  of  the  spray  nozzle  were  adjusted  such  that 
most  of  the  spray  water  entered  the  simulated  nacelle  with  the  tunnel  running. 
The  droplet -air  mixture  entering  the  simulated  nacelle  passed  through  a 
separator  box  which  contained  a  series  of  baffle  plates.  The  baffle  plates 
subjected  the  flow  to  a  series  of  sharp  turns,  thereby  separating  the  droplets 
from  the  airstream  and  collecting  them  in  a  water  film  on  the  plate  surfaces . 
Most  of  the  water  entering  the  nacelle  in  the  form  of  droplets  was  collected  at 
the  bottom  of  the  separator  box.  The  air  was  evacuated  from  the  separator  box 
by  means  of  the  suction  provided  by  the  inlet  of  a  secondary  blower.  The  air 
from  the  exit  side  of  the  secondary  blower  was  dumped  back  into  the  tunnel, 
downstream  of  the  nacelle  location.  Calibration  checks  were  made  to  determine 
the  carry  over  loss  of  water  in  the  form  of  fine  drops.  For  the  spray  flow 
rate  employed  (300  cc /sec),  the  carry  over  loss  was  about  5  percent,  with  95 
percent  of  the  water  sprayed  into  the  nacelle  being  collected  at  the  bottom  of 
the  separator.  With  a  steady  spray,  the  total  mass  of  water  collected,  divided 
by  the  time  of  collection  was  close  to  the  time  averaged  mass  flow  rate  of 
water  entering  the  nacelle  in  the  droplet-air  mixture. 

The  optical  set  up  is  shown  in  figure  1,  The  laser  sheet  of  9  mm  thickness  was 
produced  in  front  of  the  nacelle  inlet  plane  approximately  2  mm  from  the  inlet 
plane  to  avoid  illuminating  the  water  film  on  the  rounded  nacelle  entrance. 
The  illuminated  droplet  field  was  viewed  through  the  side  of  the  tunnel.  The 
camera  axis  formed  an  angle  of  approximately  20  degrees  with  the  nacelle  axis. 
In  viewing  the  droplet  field  through  the  plexiglas  wall  of  the  tunnel  at  such  a 
shallow  angle,  multiple  reflections  of  illuminated  drops  were  encountered 
within  the  plexiglas  wall.  To  alleviate  this  problem,  a  window  extension  was 
mounted  on  the  tunnel  side  wall  at  a  20  degree  angle,  such  that  the  camera  line 
of  sight  was  normal  to  the  window  at  the  end  of  the  extension. 


The  objective  of  the  experiments  was  to  compare  the  instantaneous  water  flow 
rate  into  the  nacelle  as  determined  by  the  present  non- intrusive  optical 
technique  with  the  known  time  averaged  water  flow  rate.  Ten  photographs  were 
taken  for  a  fixed  spray  rate  and  airspeed  to  yield  a  set  of  ten  successive  LWC 
measurements  in  a  steady  spray.  The  variation  in  the  drop  velocities  within  a 
given  double  exposed  photograph  or  between  two  photographs  at  fixed  tunnel 
speed  and  water  injection  pressure  was  found  to  be  less  than  5  percent. 
Therefore,  a  single  drop  velocity  was  used  for  water  mass  flow  rate 
determination.  A  sample  photograph  of  droplet  field  is  shown  in  figure  6.  The 
nacelle  boundary  was  photographed  separately  and  superimposed  onto  the  droplet 
photograph  in  figure  6.  Notice  that  the  drops  within  the  light  sheet  are 
uniformly  illuminated  and  appear  with  sharp  boundaries.  Unfortunately,  due  to 
light  scattered  from  drops  contained  in  the  light  sheet  and  that  from  the  beam 
dump  and  the  transmitting  side  plexiglas  window,  some  of  the  drops  outside  the 
laser  sheet  appear  with  faint  images.  During  image  processing  of  the  photo¬ 
graphic  negatives,  these  faint  drops,  which  are  present  outside  the  light  sheet 
may  easily  be  eliminated  by  setting  proper  threshold  criteria. 


IMAGE  PROCESSING  SYSTEM  DESCRIPTION 

Processing  system  architecture  is  depicted  in  figure  7.  Image  acquisition, 
display  and  processing  was  accomplished  using  a  De-Anza  ID-5400  image 
processing  system.  The  hardware  package  incorporates  a  vidicon  and  power 
supply  for  analog  image  formation,  three  image  refresh  random  access  memory 
channels,  RAM,  digital  video  array  processor,  and  color  video  display.  The 
analogue  signal  from  the  video  camera  can  be  digitized  by  an  A/D  converter  and 
fed  directly  to  the  array  processor  which  in  turn  controls  the  data  flow  and 
writes  the  data  into  one  of  the  memory  planes  at  a  rate  of  30  frames/sec.  The 
digitization  process  converts  each  picture  into  512*512  matrix  element 
(pixels) .  Each  pixel  is  one  byte  number  (256  resolution  level)  representing 
the  average  optical  density  in  an  elementary  cell,  the  size  of  which  dictates 
the  spatial  resolution  of  the  system.  While  digitization  can  proceed  at  video 
rates  of  30  frames  per  second,  a  program  is  used  which  creates  one  digital 
frame  from  the  average  of  64  consecutive  digitized  frames.  Thus,  the  image 
formed  has  a  low  level  of  random  noise  caused  by  the  vidicon  and  the  digitizer 
electronics.  The  digitized  image  information  is  then  stored  on  a  mass-storage 
device  for  further  off-line  processing. 

Software  residing  in  the  host  computer  (PDP  11/34)  operates  through  a  direct 
memory  access  (DMA),  interface  through  which  the  PDP-11  sends  and  receives 
information  from  the  video  processor  registers  or  from  the  RAM  channels  via  a 
driver  program.  The  vidicon  image  digitization,  averaging  and  storage 
capability  is  part  of  this  (DMA)  interface  software. 

To  study  a  droplet  picture,  the  negative  is  taped  on  to  a  transparency 
containing  a  1  cm  x  1  cm  grid  formed  by  fine  lines.  The  grid  lines  are 
carefully  oriented  with  respect  to  the  droplet  images  on  the  negative  such  that 
none  of  the  drop  images  is  intersected  by  the  grid  lines.  The  negative, 
together  with  the  grid  lines,  is  then  mounted  on  a  flat  light-table.  The 
vidicon  is  focused  on  the  plane  of  the  background  illuminated  negative  by  means 
of  a  macro  lens.  The  magnification  on  to  the  vidicon  is  adjusted  such  that  a 
1  cm^  area  of  the  negative  enclosed  by  the  superposed  grid  lines  nearly  fills  a 
video  frame  containing  the  512  x  512  pixel  array.  After  accounting  for  the 
photographic  magnification  of  the  image  on  the  negative  and  a  specification  of 
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4  pixels  for  the  minimum  size  drop  accepted,  an  overall  resolution  of  0.1  mm 
was  estimated  for  the  diameter  of  the  smallest  drop  measured. 


The  negative  was  scanned,  one  grid  square  at  a  time  over  the  projected  frontal 
area  of  the  nacelle.  The  calculated  drop  volumes  were  then  added  up  for  all 
the  squares  to  evaluate  the  liquid  water  content  over  the  projected  frontal 
area  of  the  nacelle.  This  procedure  ensured  that  every  drop  is  accounted  for 
and  is  counted  only  once.  An  example  of  a  digital  subimage  representing  one 
grid  square  is  shown  in  figure  8.  An  enhanced  image  with  defined  droplet 
boundaries  is  shown  in  figure  9. 

Image  processing  software  was  developed  to  detect  drop  edges,  define  droplet 
boundaries,  and  calculate  the  area  of  drop  images.  The  excellent  contrast 
between  the  images  of  drops  within  the  light  sheet  and  the  background  allowed  a 
simple  thresholding  criterion  to  define  regions  containing  drop  images. 

The  various  image  processing  programs  along  with  computer  program  needed  for 
droplet  statistics  determination  are  attached  under  Appendices  A  and  Bj 
respectively. 


RESULTS  AND  DISCUSSION 

Results  are  presented  here  for  fixed  values  of  spray  flow  rate  and  airspeed  in 
the  wheel  spray  simulation  tunnel.  The  flow  rate  through  the  spray  nozzle  was 
maintained  constant  at  0.3  liters/sec,  which  resulted  in  a  water  jet  velocity 
of  22  m/s  based  on  the  total  flow  area  of  the  seven  tubes  in  the  spray  nozzle. 
As  discussed  later,  this  flow  rate  was  chosen  to  simulate  extreme  cases  of 
engine  water  ingestion  resulting  from  wheel  spray.  The  airspeed  was  maintained 
at  61  m/s.  The  large  difference  between  the  speeds  of  the  air  and  the  water 
jets  caused  an  aerodynamic  breakup  of  the  water  jets  into  small  drops,  which 
accelerated  along  the  flow  direction  to  approach  the  airspeed.  The  measurement 
station  was  located  3  m  downstream  of  the  spray  nozzle  and  immediately  upstream 
of  the  simulated  nacelle  inlet. 

The  droplet  velocities  were  determined  from  a  double  exposed  photograph  of  the 
moving  droplet  field  as  shown  in  figure  4.  A  100  ns  time  interval  between  the 
two  pulses  was  employed.  The  droplet  velocities  were  found  to  be  within  3 
percent  for  a  large  number  of  droplet  pairs.  Furthermore,  within  this  small 
variation  of  droplet  velocities,  no  correlation  was  found  between  the  drop 
velocity  and  size.  A  single  average  value  of  velocity  was  assigned  to  all 
droplets.  This  average  value  was  found  to  be  only  41  m/s,  i.e.,  67  percent  of 
the  airspeed.  Thus,  the  3  m  distance  between  the  spray  nozzle  and  the 
measurement  station  was  insufficient  to  accelerate  the  drops  to  the  tunnel 
airspeed . 

A  series  of  ten  photographs  was  taken  during  the  period  January-February  1986 
for  the  liquid  water  content  determination.  The  width  of  the  light  sheet  was 
maintained  at  9  mm.  The  photographic  negatives  were  analyzed  by  the  image 
processing  technique  discussed  in  this  report.  The  volume  of  a  drop  was 
calculated  as  the  volume  of  a  spherical  drop  having  the  same  equivalent 
diameter.  The  equivalent  diameter  was  calculated  from  the  image  area  of  each 
drop.  For  a  non-spherical  drop  this  procedure  leads  to  a  higher  value  of  the 
calculated  drop  volume.  The  percentage  error  introduced  therefore  depends  upon 
the  how  different  the  drop  is  as  compared  to  the  spherical  shape. 
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The  results  of  the  image  processing  analysis  are  shown  in  table  I.  The  raw 
droplet  data  has  been  shown  under  Appendix  C.  For  each  of  the  ten  droplet 
pictures  analyzed,  the  total  number  of  drops  counted  over  the  nacelle  frontal 
area,  the  mean  droplet  diameter,  the  volume  weighted  mean  diameter  and  the 
total  liquid  volume  are  presented.  The  standard  deviation  of  the  mean  and 
volume  weighted  mean  droplet  diameters  as  calculated  from  the  ten  pictures  was 
in  the  range  of  7  to  8  percent  of  the  average.  This  indicated  that  the 
picture- to-picture  variation  of  the  calculated  mean  diameter  was  ^^.latively 
small.  The  pic ture - to -pic ture  variation  of  the  total  liquid  volume  was  larger: 
the  standard  deviation  of  the  ten  measurements  was  31  percent  of  the  average. 
The  average  value  of  the  total  liquid  volume  from  the  ten  pictures  was  used  to 
calculate  the  average  liquid  water  content,  i.e., 

(LWC)a  -  '°w^w  —  450  gm/m^ ;  where  ^w  -  1.0  gm/cm^ 

Af  .& 

The  average  mass  flow  rate  of  water  v*as  then  calculated  as 

™w  -  ( LWC ) a  .  Af  .  Vd 
-  332.5  gm/sec 

The  average  mass  flow  rate  of  water  in  the  moving  droplet-air  mixture  as 
measured  by  the  present  non- intrusive  optical  technique  thus  indicates  a  flow 
rate  approximately  10  percent  higher  than  the  actual  value  of  0.3  1/s.  Reasons 
for  this  discrepancy  between  the  measured  and  the  actual  mass  flow  rates  are 
discussed  below.  It  should  be  noted  however,  that  there  is  a  significant 
picture- to-picture  variation  in  the  instantaneous  water  mass  flow  rate 
determinations  caused  by  measured  variations  in  the  instantaneous  LWC  values. 
The  test  droplet  spray  was  quite  dense  in  comparison  with  the  certification 
requirement  of  4  percent  by  weight  of  liquid  water  in  the  airstream.  The 
present  0.3  1/s  water  ingestion  rate  in  a  61  m/s  airstream  over  the  frontal 
area  of  the  15.24  cm  diameter  nacelle  translates  to  a  water  flow  rate/airflow 
rate  ratio  of  22.5  percent.  Such  a  high  water  spray  rate  was  used  to  simulate 
extreme  cases  of  engine  water  ingestion  resulting  from  wheel  spray. 

There  are  three  factors  that  affect  the  accuracy  of  LWC  measurement  by  the 
present  technique: 

1)  Droplets  largely  outside  the  light  sheet  but  at  the  boundary  are 
partially  grazed  by  the  light  sheet  and  show  up  on  photographs,  thus 
increasing  the  LWC  measured. 

2)  The  procedure  for  calculation  of  drop  volume  from  its  non- spherical 
images  on  the  photograph  relies  on  an  equivalent  diameter  which  is 
calculated  from  the  enclosed  area  of  the  image.  This  procedure  tends 
to  over-predict  the  volume  of  the  non- spherical  drop  and  hence  leads 
to  a  higher  value  of  measured  LWC. 

3)  The  illuminated  drops  contained  in  the  light  sheet  are  viewed  by  the 
camera  through  a  dense  spray.  Therefore,  there  is  a  possibility  of 
some  illuminated  drops  be  .ng  masked  by  droplets  present  in  the  view 
path  of  the  camera.  This  masking  will  result  in  a  lower  value  of  the 
measured  LWC. 


TABLE  I 

Summary  of  Data  Photographs 


Va  “  61  m/s; 

;  Vw  -  21.7  m/s; 

Qv  -  300 

0 

cm  /sec ; 

x  - 

3  m; 

Af  —  182.4  cm^ 

S  =  9  mm 

Negative 

No.  of  Drops 

Mean  Dia. 

Vol. 

Mean 

Dia . 

Total  Liquid 

No . 

- 

mm 

mm 

3 

mmJ 

1 

267 

0.55 

0.89 

97.4 

2 

207 

0.58 

0.86 

68.2 

3 

256 

0.51 

0.87 

89.0 

4 

155 

0.58 

0.86 

52.3 

5 

179 

0.63 

0.92 

88.0 

6 

212 

0.56 

0.88 

75.0 

7 

243 

0.59 

0.91 

95.7 

8 

128 

0.50 

0.68 

21.5 

9 

337 

0.53 

0.80 

91.8 

10 

169 

0.60 

0.87 

57.8 

Average : 

0.56 

0.85 

73.7 

Standard  Deviation: 

0.04 

0.066 

23.0 

Std 

Dev.  percent  of  Mean: 

7.0 

7.8 

31.0 

NOTE : 

„  n.  Zn:Di 

Mean  Dia.  - 

;  Vol.  Mean  Dia. 

ZniDi4 

Zn^ 

Zn^3 

Total  Liquid  Vol.  -  Zrw  n  D,-3 


Despite  these  three  sources  of  errors,  the  present  non- intrusive  technique  is 
still  considered  a  good  method  for  measurement  of  the  liquid  water  flow  rate  in 
a  moving  droplet  -  laden  airstream. 

The  droplet  size  distribution  determined  from  the  more  than  2000  drops  counted 
in  the  ten  pictures  is  shown  in  figure  10. 


CONCLUDING  REMARKS 

1)  A  non  - intrus ive  optical  technique  has  been  developed  for  quantitative 
determination  of  instantaneous  liquid  water  mass  flow  rate  in  a  droplet 
laden  airstream.  The  technique  is  generally  applicable  to  the  problem  of 
water  ingestion  into  an  engine  resulting  in  up  to  22%  water  by  weight 
(rain  or  wheel  spray). 

2)  The  technique  yields  instantaneous  spatial  distribution  of  droplets  at  the 
nacelle  inlet  plane  as  well  as  the  droplet  size  distribution. 

3)  Significant  variation  in  the  instantaneous  values  of  the  liquid  water 
content  was  encountered  in  a  droplet  laden  airstream  produced  by  injection 
of  a  steady  water  spray  in  a  steady  airstream.  The  standard  deviation 
among  ten  separate  instantaneous  determinations  of  the  LWC  was  31  percent 
of  the  average  LWC. 

-0  The  average  liquid  water  flow  rate  as  determined  from  the  average  LWC  and 
drop  velocity  measurements  was  approximately  10  percent  higher  than  the 
actual  spray  flow  rate. 
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ACCEPT  101, ANSA 
IP  (AHSA.EU.  ' Y '  1  GO  TO  950 


0 

C  F  I  rat  detector 

i; 

1 1 1 


0044 


TYPE  *,'  Ucgtn  coarse  detection' 


FORTRAN  IV-P  J  Vn2-r.l  10:00:17  2G-APR-0G  I  2 

DROPS. FTN  /  riCOLOCKS/VR 


no. ir. 

CAM.  SETUP 

»)iM  •> 

TYPE  Smoothing  before  compute  gradient  <Y/N)' 

nm  / 

ACC  LPT  mi. ANSA 

nmii 

IK  < ANSA. EO. ' M ‘  )  CO  TU  4 

Wlf) 

CALI.  SHOOTII 

nnsn 

so  ro  ti 

DOS  1 

4 

1)0  1  K-  1 , 2SG 

r»0'i.r 

J 

t‘/I.T<  K  )«K/2 

DOS.) 

|  * 

CALL  VL TCUIU  1  .  I VLT  ) 

i: 

I* 

Save  original  picture 

nos  4 

•i 

CALI.  SC 7 UP 

DO'j'j 

OIM:h  ( UN  If" /.TYPE- 'SCRATCH ' .ACCESS- 'D IRECT ' . INITIALS I7E 
i  records  izc  - 1  ;:u > 

nn'»r, 

NULL-  1 

1)00  / 

DO  in  K-1.J2 

noon 

i.  in  ■<  k-  i  )*ni+ 1 

0000 

ICOII-  I  MACIi  <  n,  1.  IN,  1G,  'R'  .PIXEL  I 

POOD 

I SCAT -J CSC  <  7.1*  IXEL.NDLK.4nsG.  1  ,  NBLK  ) 

000  1 

in 

i; 

i; 

C 

CONTINUE 

Transfer  smooth  picture  to  channel  n 

nor,.? 

CAM.  SETUP 

000  ) 

call  !Hsi:i.<2.n,n) 

001)1 

CALL  OUTNtll. <  1  ) 

0))00 

C 

i: 

c 

CALL  KHOPER 

Compute  gradient 

nor.o 

CALL  SETUP 

000  / 

c. 

0 

c 

CALL  CHAD  IE 

Re-Jtore  original  picture 

noon 

CALL  SETUP 

0000 

N11LII-*  1 

oo  fa 

do  ;;•)  K"  i .  32 

00/  1 

1  I N-i  1C-  1  )  1 

oo/;? 

I  ST A7-.IL EC <  /  .  P  I XEL  , HULK  ,  4090  ,  1  .  NQLK  ) 

00  /  t 

ICON- I  HAUL  <  n . L IN, 1 G .  'V' .PIXEL) 

00  M 

zu 

i* 

CONTINUE 

i; 

c 

Threshold  gradient  picture 

rT07r> 

CALL  SETUP 

00  /l> 

CALL  IHSEU  1,0.0) 

00  /  / 

CALL  OU 7 NHL <  4 ) 

00  / 0 

CALL  KHOPER 

00  /•) 

CALL  SETUP 

0000 

so 

TYPE  ■* .  ’  Threshold  gradient  Image  '.  ICORTA 

non  i 

CALL  SETUP 

non.; 

DO  SI  JET™ 1 , ICORTA 

non  i 

s  l 

ivl  n  Ji.  r  )-n 

000  1 

DO  'G2  JLT-  ICORTA*  1 ,25G 

1/7:00?  17 


2G-APR-05 


FORtJUlf  IV-  S  V/T2-0I 
DROPS. FTM  /  fit  :UL0CK3/WR 


nuns 
00. ‘Hi 
0(111  / 
mum 

ilDII'l 

ODD  t 
00*1  ! 
Illl'l  I 

on1)  I 

DifCi 
oner, 
•»>»•»  > 
U 


oom) 
o  10.0 
0  10  1 
oio.! 
o  I  0  t 
0  10  1 


01  O’i 
0100 
o  1  o  / 
0  100 
0100 
0  1  to 
0  111 
on;! 
out 
out 
o  1 1  •> 

Olh. 
0  11/ 
•mu 
ouo 


o  t  r >o 
0  12  1 
o  i ;!;? 
o  i ,:i 
o  i  ;•  1 
0  1  »!’> 


n  i r, 
o  i > 
0 1  jn 


101 

(jJJ 

70 


C 

i: 

c 


no 

i; 

i; 

C 


or, 


i.?o 

t: 

C 

r; 

mi 


I VI.  T<.1LT  l»20r. 

CAM.  VIDTIIR1  1  .  1 VI. T  > 

TYPE  Accept  throshold  <  y / n  )  * 

ACCEPT  101. ANSA 

IE  (ANSA. ME .Ml  CO  ft)  GO 

TYPE  * .  '  Enter  gradient  Imago  threshold’ 

ACCEPT  *.ICUUTA 

co  ro  'Hi 

I"  OR  HAT  (  A  1  1 
CAM.  SETUP 

DO  70  JET  »  ICORTA-*  1 , 2GG 

i  vi.  r<  jl  r  i  =.ii.  r- 1 

CAM.  VI.  rcoim  .  I  VL  T  ) 

CALL  SETUP 


Save  gradient  picture 
HULK*! 

DO  110  K ®  l  .32 
I.  IN  «<  i:-  1  1*  Hi*  1 

I  COI)-»  IHACE  <  I  .1.  IN.  16.  R  1  .PIXEL  > 

IS  f A  f -.IE SC <  /  ,P  IXEL  .  MBLK  ,  HU'Jlj .  1  .HULK  I 
CON  f  INUL 

Humovo  background  variations 

TYPE  Remove  background  variations  !Y/N> 

ACCEPT  101. ANSI) 

IE  <  AIISO.  HE.  'll'  )  Cl)  TO  D!i 
CAM.  I  MS  El.  <2.0.0) 

CAM.  OUT  Hill.  <  4  ) 

CAM.  I  HOPEH 
CAM.  SETUP 
CAM.  INTEL  <  1 .0,0) 

CAM.  OU7HIIL<2) 

CAM.  E  HOP  EH 
CALI.  SE  TUP 

co  rn  inj 

1  SI)-*  1.00  I  Variations  scale 
CAM.  II ICPAS1  ISO,  1  ,0  ) 

CALL  SETUP 

Kostore  gradient  picture 
NI1LIM1 

1)0  1 20  K»1 .12 
L  1 II'* v  K -  I  )*1E,*1 

IS  f AT « JLEE 1 /.P 1XEL . HULK, 4006. I .HULK  ) 

I  COD’*  IHACE  I  2  .  L  1 H  .  1  G  ,  'V  .PIXEL  ) 

CONTI  HUE 


Compute  Histogram 
TRim-O 

fYPE  * . '  Computing  mod  If  tod  PDF' 
DO  'JO  »:='0,2’j'J 


iy-  S  V03-51  10:50:1?  2G-APR-05  € 

DROPS. FTM  / fU :UL0CKS/WR 


012*1 

■JH 

ii  is  n  i:  i’0. 

/T  1  Ilf 

Ill)  150  K-  10.500 

Hill 

icoim  ihai;i:  <  i .  i: .  i . '  R  1  .hup  ic  ) 

HI  12 

tcuii-  imacei  2.i: .  l .  •«  ’  .uucka  i 

m  n 

110  150  .1 "  10, 500 

HIM 

II  <  I1U0RAI  3)  .IU.2CR0  )  00  TO  150 

HI  I’i 

A"  I  V<  IIUCRAI  J  I  I 

H  1  tii 

I  I  ™  1 VI  IliJI*  1C<  J  I  I 

HI  t  ' 

ii i ■;  r <  1 1  > -n is  r 1 1 1  »♦  i  i«a 

0  1  JIJ 

150 

I  * 

CUN  f  IHUL 

1  # 

<; 

c 

Display  hlsto<jran  and  threshold  determination 

uni 

UMAXO-oJ. 

it  i  in 

Hi)  1  HH  K»0. 255 

HIM 

11  (IHirtKI.OT.OHAXO)  DHAXO-H  I3T< K > 

HI  <2 

100 

con  r  ihiji: 

Hill 

III)  IIJ.2  K '*0.255 

HIM 

ii c; rt  KHiism  i/ohaxo 

H  1  15 

1U2 

continue 

H  l  Hi 

CAI.L  1 N5FI.  <  1  .  H  ,H  I 

H  l  1  / 

CAI.L  OH  f  Hill.  <  4  ) 

H  Mil 

CAL  1.  FMUPER 

HI  10 

CALI.  5L  riJI* 

H  1  50 
HIM 

H  I'i.’ 

H  1 5  1 

1 50 

1)0  ion  ic-|.25G 
ivl  rt  ic  )=h 

CAM  V  ID TIIRI  1  .  IVLT  ) 

I  11111=0 

HIM 

HI'i'i 

CAM.  I)  1 511 1 3<  II  1ST  ,  ITHR  ) 

r  VI*  1!  ■* ,  1  Enter  (1)  simple.  (2)  blmodal  distribution’ 

H  I'Hi 

ACCEPT  * ,  I MODE  5 

H  1  5 

CAM.  3EUISO<  HIST.  ITIIR.  I MODE 5  ) 

H  1 50 

300 

TYPE  Estimated  tliraaho  Id  ITIIR 

H  1  5  0 

CAM  SETUP 

H  Hill 

1)0  10!  JET"  1 .  ITIIR 

H  I'i  1 

301 

IVL  M  .IL  r  >"0 

H  Hi.* 

1)0  IH.*  JET"  ITIIH+1 .25G 

Him 

302 

ivl  rt  .ii.  r  1-255 

Dili  1 

CALL  VI.  rcORIl  .  IVLT  ) 

H  lli'i 

call  si:  rui* 

H  1 •>  >i 

i  vi.  r<  i  i  mj 

H  1  5  ' 

IVI.  f  {  255  1-  l 

H  1  liH 

CALL  VI.  fCUK<  1  .  IVLT  ) 

H  Hi’l 

cam.  si:  nil* 

H  1  /H 

CALL  IIOUHEAt  IVL  f  ) 

HIM 

CALL  SETUP 

H  i  »:* 

cam.  nismstiiisr.  ithr) 

HIM 

TYPE  Accept  threshold  TY/MI' 

HIM 

ACCEPT  101. ANSA 

H  l  *5 

II  1  ANSA  .  HE  .'H  I  CO  TO  350 

H  1  .Mi 

TYPE  Enter  threshold’ 

II  l  f  ’ 

ACCEPT  *  ,  I  TIIK 

H  1  Ml 

11  1  ANSI) .  HE  .  'll’  )  CO  TO  320 

H  1  M 

CAM.  I  MSI  1.  T  1  ,:l.U> 

H  l  HH 

CAM.  D'lfNIll.  (2) 

HIHI 

CAI.L  IIIUPLR 

rOKTIWN  IV- 
DHOPS.FTM 


S  v/T2-r.  i 

/ f  K : UL0CK3/VR 


1(1:00:  17 


2G-*rn-o5 


CALL  Si:  TUI* 

tin  ro  jiw 

C/M. L  IIHiPASI  ISO.  1  .01 

cali.  si:  mi* 
co  ro  jiw 

r Y l» Default  region  tho  entire  picture:  change  (V/NI 
ACCEPT  IIIl.AHSA 
IK  IAMSA.L0.  'VI  tiO  TO  400 


HI.  ■■■'>  1  0 
n  ■;=•'>  ii) 

«;*i  ro  400 

r VI’i;  Die  cursor  to  define  region/push  white  key 

CALL  KLUIOCJL.SS.HL.NS) 

AH  I  HAH -ARM  I H 

TVIM:  Hog  Ion  definition  SL-*.3L,'  S3-', SO,'  ML3'  , 

I  vi.  r<  i  i-o 
i  vi.  r< ,? 

I VI.  TCVift  l-.l'iO 
CALL  VI.  rC0K(  1  ,  1VL  r  ) 

CALL  :;l  rui* 

r V l» i:  Computing  dropt  size' 

CAM.  (in  IASI  1  ,SL  .  JS.HL .NS.AUUK  .POUF  .LIMIT.  NUMODJ  ) 


call  :;l  ru 
tv i»i: 

CAM.  (ill  r A 

i  vi.  ri  1 1  mi 
i  vi.  r<  :*'ic,  i 
cam.  vi  rc 

CAM.  H'llUl 

cam.  -;l  ru 
tvpi:  *.• 


vi.  rcoR<  l .  i  vi.  r  i 
IlMHIILA'  IVLT  I 
■;l  rui* 

Lnd  coarac  detector' 


CLOSHUHl  f-/  ) 

Save  ro3ult3  on  disk 
HUHDIMiT 

■;cai.l::33calc*3cai.k 

rvpi:  DATAIIA  ii;:  Initial  element  ’  .KRC DC 

no  Min  K=*  1  ,  HUMOUJ 
HUHIlIt  t  HUMIIU  *  1 

t:iu:i)c->i;ui:iic  »i 
PI. “L  IMI  rt  I  .K  * 
pc»i.  mm  ;:.k  i 
Cl. "1.  Ill i  r<  :i.k  i 
cc =i.  in  i  r  f  4  .k  » 

AIILA-AIIUI  1 1C  >* SCALE 2 
IMH  IM-*P1IUF(  K  I* SCALE 

i/iu  n.  r  I'lcm.wci  I  rukE 

II  (  ICKLDC.EU.  12)11  I  00  TO  2025 
com  r  imiji: 

rvpi:  DATAI1A3E  :  Final  e  lement ’ , KRF DC- 1 

TVIM;  *,•  Hunlier  of  ol)Joct3  In  thli  f  mine  '  .  NUMDR 

lilt  1  fl'.l  T  '  1  I  1 IHIFU  I  Update  general  Information  record 

Cl MSlt UHIT-J I 

Mil) 


’  V  V 1 


FORTRAN  * IV-PL  V82-* 51  14)89:82  25-APK-85  PA  1 

COTAS. FTN  /TR : OLOCKS/VR 


nnm 

c 

c 

c 

min?. 
non » 
non  ( 
nuns 
nnns 
nnn ' 
nnnn 
nnnn 
im  in 
nn  l  i 

i)i)  12 

C 

c 

c 

nn  1 1 
nn  i  i 
nn  is 
nn  i  ■> 
nn  i  / 
nn  i  n 
nn  i  •» 
nn.sf 
Du?  i 
nn,’,: 
nn,:  l 
nn  ■ » 
nn.’S 
H’l.’ii 

nn.: ' 
nn.:  n 
nn.”> 
nnjn 

C 

i; 

c 


SUORUUTINE  GOTASICHAN.SL.SS.NL . NS . AUUF . PBUF .LIMIT .NUMOBJ  ) 
Drops  goonetrlcal  caractar 1 *ot I  on 

mrrcrR‘2  chah.sl.ss  ml. ns 
IIY  rr,  P  1 XC L ( 5 1 2  .  Hi) 

COMMON  /1'IMTA/PIXI !L 
I N  ft:  lit:  11*2  FLAG!  S  1 2 . 2  ) .  ID<3lT) 

IN  ri:r;ru«2  ISLSI  Wl .  (STS!  9/1)  .EL<  98)  .FS<  ’ill) 

in ti:i;i:r«2  cm th (  nn )  ,uluin< *jd  .2 ) .  ilmik  in.  2  ; 
u.  rr<;i:u*2  DLUtn.iirviD 
ihilglum  hpixscjd) 
real  *4  npritcM) 

IM  AI  M  Aliur  <  1  I  .l‘liur<  1  ) 

IN rtGLR*2  L I H  IT  <  4 .  1  > 

Initial  1  za 

CALL  2  I A  (  ri.AG<  1,1).  1824  ) 

CALL  2  1A<  IM’l  ll  .  inn  I 
CALL  2  I A<  HI*  1  XT  .  10//) 

CALL  XIAILL.nifl 
CALL  2 lAt LS .98) 

CALL  2  lAfCHTUf  1,1),  l  IDT  I 
CALL  2  I A  <  II  i:  •;  1  N<  l  .  1  >  .  HOT  ) 

CALL  2  IAI  U  NIX  l  ,  1  >  .  1UD  ) 
call  i  r  ia<  i2nnn,  isi.:;, no  i 

CALI  1  r  1 A <  J2DDO.  IG3U.9U  ) 

H'JI4URJ"8 
i  sw  <  l 
I4LWC*8 
HOI  on 
01  I)  II)’-  1 

ms  -ns 

IISL- SI-1 
MLL-HSL+ML 

Main  loop 


nm  i 

i»o  tin  i.-i  .ml  .  tr» 

nn  i:: 

1  lll.OCKaM  1 N04  10  ,  IIL-L  > 

nn  n 

1.  IIH>1  >H5I. 

nn  M 

1  CUD  -  I MAUL ( CHAN , L INC . IOLOCK, ’R ' .PIXEL ) 

nn  is 

DO  HID  ICK- 1.10 

nn  is 

H  1  "SS 

nn  i  / 

1.  II4IA-L  INC  ‘KK-1 

nn  Ml 

MIWC»W  1  Start  searching  for  an  object 

nn  in 

185 

DO  llll  .l-Ml.NS'i  1  look  for  left  od’jo 

nn  in 

ii  (  p  i :<i:l f  j.ki:  )  i  go  ru  128 

DIM  1 

i  in 

cum  r  inui; 

c 

Mo  worn  cd'joa  on  tills  lino 

nn  42 

<;o  ru  im 

nn  1  1 

i2n 

II  1  ->.)  1  I  nf  t  edge 

nn  i  i 

Do  IJ'f  .1-141, MSG  I  Look  rtgh  edge 

i)i)  IS 

II  !  .HOf .PIXLLI J.KKI )  GO  TO  14 8 

nn  in 

13  8 

CUH  r  1 HIJC 

-A6- 


u  -.no az 


25-APR-00 


pa  e 


fOttMN  IV-PL  V02-31 
COTAS. FTN  TIUULOCKS/WR 


not  ’ 

Amu  no 

c 

c 

c 

ooio 

Ild'iO 

nos  i 

c 

c 

ruin-' 

i)  i)  'j  4  10! 

c 

c 

0030 

c 


OOOG 


003 1  1  (12 

OO'iH 

DOS') 

OiH.O 
000  t 


0002 
000  I 
000  4 
ooo'i 
OO'iO 
iino  / 


ooon 

IllUlJ 


i; 

i; 


J-1433*  1  I  High  edge  13  the  picture  boundary 
I42-J-  l 

Find  10  of  this  object  from  provloui  flag  lino 

CALL  HA  fCHt  FLAG!  1,3-1  3W  )  ,  HI  .  H2  .  HUM  ,  I U  .  1433  ) 

HI  win- !0<  1  ) 

IF  UIUH.LT.  2)  GO  TO  101 

If  innv  II)  condense  them  Into  one 
CALI  COM  I  Of  ItPLIt  .HP  1X3  .  I'JU  .  1 33  3.  EL  .  C  3  ,NUH  .ID) 

on  rn  i#2 

If  4  HUN. HI. 0)  GO  TO  102 

If  thrern  In  not  ID  them  creates  one 
CALL  I  1  HIM  90.MLW1D.MP  IaS.ULDU) ) 

bo.  jin  count  lug  per  I  note  r  now  object 
It  I*  L  It  (  ML  W I U  I-I42-M  1  +2 


And  update  flag  buffer  and  counters 

CAI  I.  I  r  I  At  HLWil).l  LAC<  Ml  .  I  GW  I .  H2-N  1  *  1  1 

Ml  WC’MLWC*  l 

cm  f it (  nr wc  .  isv>*Mi:vto 

lllOIMt  ml  wc ,  ISUK-Ml 

ItMDt  MLWC  .  !  jW  )»I42 


Update  statistics  buffer 

Ml*  1X3!  MLWIO  >“KP  1  :<!»<  ML'W  1 1)  1  ♦  1 

II  <  l3I.S<MFWIH>.i;r.UMIA»  131.34  nlwid>-l  I  m  I  a 

II  4  13334  MEWIU >  GT .Ml  )  IS3S<  MCWIin-Ml 

IK  1 1  34  MGWII)  >  .4. T  .M2  >  ESI  NLWIl)  >*M2 

LI  (MLWID)-l  1 14 1 A 

CALI  I'l  It  IN4  MULC.Mr.WC.CH TR4  1  .  I OW  > .  HGV1 1>  .  UPGR  ,  HUG  I  MI  1 ,3-ISV), 
1  BEG  1144  I . ISW). ICN04  1 ,3- ISW). ILM04  1 . 13WI  ) 


'lend  back  for  more  data  on  the  same  line 
M 1  -142  >2 

II  < Ml. IE.. *33)  GO  TO  103 


Search  buffer*  for  terminated  drops 


00 70  103 


00  . » t 
00/2 
Oi*  >  I 

"O' I  100 
■to 


CALL  LI40I  I.  14  ItPKIt  .NP  ITS.  1SLS.  1033. El.  ,  E3  .  MI.L  .  L  t  N  I A  .  C14TR  4  1  .3-I3W). 

1  CM  f K  <  1 .  I  >W I . MOLC , MLWC . MUMOUd  DKG  I  MI  1 . 3- I 3W I , BEG  I M4  1 . I 5W I 

2  1 1 NU<  1  .3 -ISW)  .  | UNIX  1  .  ISW)  .All  UK  .IMIUK  .LIMIT) 

CALL  ZIA4 FLAG!  1 .3- ISW)  .312  )  I  Zero  flag  line 
ISW"3-I3W  I  Switch  the  flag 

HOLC-MI  WC 

CUHMIIUC  I  Close  wain  loop 
III  ruitH 


FORTRAN  IV-F 
HATCH . FTN 


V02-51 

/TR : BLOCKS/VR 


If i»l :J0 


04-APR-B5 


ViV 


Kirtv-CvM  - 


SUBROUTINE  HATCH ( FLAG .N 1 , N2 . NUN , I D . NS  > 

RETURN  1.0.  OF  OBJECT  FOUND  ON  PREVIOUS  LINE 
FLAG  -  OLU  FLAG  BUFFER 
N 1  -  BEGIN,  N2  •  END  OF  OBJECT 
HUM  -  NUMBER  OF  OBJECTS  LOCATED 
ID  -  I.D.  OF  OBJECTS  LOCATED 
INTEGER  FLAG! 1 ) 

IM7EGER-4  ID< 1 > 

Ml  -  HI  -  I 
M2  -  M2  *  1 
IF  (Ml  .LT. 1  )  Ml  -  1 
IF  (M2.GT.NS)  M2-NS 
MUM  -  a 


FIND  ALL  FLAGS 

00  10  J-Ml .M2 

IF  <  FLAG( J  I.EQ.K)  GO  TO  1* 
IF  ( FLAG<  J  I.EQ.0)  GO  TO  IB 
K  -  FLAG< J  ) 

NUH  -  MUM  ♦  1 
I IX  MUM  )  -  K 
1  CONTINUE 

IF  (MUM.LT.2I  RETURN 

REJECT  DUPLICATES 

CALL  CORSRT < 10, ID.NUMI 
N  -  1 
K  •  10(11 
DO  20  J-2.NUM 
IF  (K.EO. ID( J ) )  CO  TO  20 
M  -  M  I 
IIXMI  -  I  D(  J  ) 
r  COM r I  HUE 

HUM  -  N 
HI  TURN 
END 


.< -  •r.t  -T,  V  /  <rm  'r+  'r<  ^ 


FORTRAN  IV-r  V02-51  19J0H04  04-APR-BB  P  1 

PERJM.FTN  /TR : BLOCKS/WR 


0001 


0#02 

0003 

0004 


SUOROUTIME  PER !H( NOLC .NEWC ,  NECNTR , NEW ID .RPER.OBEGIN.NBEGIN, 

1  OEND.NEND  ) 
r 

C  THIS  SUBROUTINE  ACCUMULATES  PERIMETER  MEASUREMENTS  FOR  THE  'OLD' 
C  LINE  BETWEEN  N1  AND  N2  OF  THE  'NEW  LINE. 

C 

INTEGER  NECNTRI  1  >,0BEGIN< 1  ) . NBEG  INI  1  )  ,OEND<  I  )  , N  E  N  D (  I  ) 

REAL  *4  RPER<  1  ) 

IF  (NOLC.EO.0)  RETURN 


C 

C  SEARCH  OLD’  LINE  FOR  MATCHING  OBJECT  SEGMENTS. 
C 


0005 

0006 

0007 

0000 
0009 
00  10 
001  1 

0012 


LI  -  0 

DO  100  J-l  .NOLC 

IF  ( UENIX  J ) .GE .NBEG l N( NEVC  1-1 .AND. 

1  OBEG I NI J  > . LE .NEHD<  NEWC )+ 1  )  GO  TO  101 
GO  TO  100 

101  OHEWC  -  NEWC  -  1 

IF  (ONEWC.LT. 1 )  GO  TO  301 

IF  ( NEWID . EQ . NECNTRI OHEWC ) .AND .OEND<  L2 ) .GE . NBEG I N< NEWC ) ) 

1  GO  TO  302 
CO  TO  301 
C 

C  IF  PERIMETER  HAS  BEEN  ADDED  FOR  PREVIOUS  PARTICLE  SEGMENT  ON 
C  MEW  LINE  WITH  SAME  I.U.  AND  THE  OLD'  MATCHING  SEGMENT 

C  OVERLAPS  BOTH  MEW  SEGMENTS.  THEN  THE  OVERLAP  IN  THE 

C  CURRENT  SEGMENT  IS  SUBTRACTED  FROM  THE  PERIMETER  VALUE. 

C 


0013 

302 

RPERI NEWID )  -  RPER(MEWID)  -  SORT ( <  OEND( L2  > 

1  -  FLOAT  <  MENU ( ONEWC  )  >  )**2  ♦  l.»  ♦  NBEG I N( NEWC  >  - 

0014 

GO  TO  303 

0015 

301 

RPEIUHEW101  -  RPE1UHEW1D)  •*  SORT!  <  NBEG  IN<  NEWC  ) 

1  -  FLOAT ( OBEG 1 N<  J  > >  >*«2  ♦  1  .  ) 

0016 

303 

LI  -  J 

0017 

K  -  J  *  1 

0018 

GO  TO  102 

0019 

100 

CONTINUE 

0020 

102 

IF  ( L 1 . EQ .0  )  RETURN 

002  1 

L2  -  0 

0022 

IF  (K.CT.NOLC)  GO  TO  200 

0023 

UO  103  J-K.NOLC 

0024 

IF  (OENU< J 1 .GE .NBEG IN( NEWC  >-l .AND. 

1  OBEG I N( J ) . LE . NENU( NEWC  > ♦ 1  I  GO  TO  104 

0025 

GO  TO  103 

0026 

104 

RPER<  NEWID )  -  RPERI NEWID  >  *  (OIEGIR(J)  -  0END<J-1>> 

0027 

L2  -  J 

0028 

101 

CONTINUE 

0029 

200 

IF  (L2.EQ.0)  12  ■  LI 

0030 

RPERI NEW  10 )  -  RPERI NEWID  >  ♦  SORT! ( NE NO < NEWC ) 

1  -  FLOATIOENUI L2 ) >  »**2  ♦  I.* 

RETURN 

END 


NEND<  ONEWC > 


003  1 
003  2 


18159)52 
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Fortran- iv-  s  v02-5i 

CON ID.FTN  /TR  < 8L0CKS/VR 


.  1 


0001 


0002 
0003 
0004 
0000 
000G 
0007 
0008 
0009 
00  10 
00  1  1 
00  1  2 
0013 
00  1  i 
0015 
00  1G 
00  17 
00  11) 
00  19 
0020 
002  1 


SUBROUTINE  CON ID< RPER . NP IXS . ISIS . ISSS .EL ,ES .NUN . ID > 

C 

C  USED  TO  CONCAT INATE  MANY  I.O.'S 
C 

INTEGER  OLD  ID. ISLS< l  >. ISSSI 1),EL<  1  I.ESt  1  ) 

INTEGER**  NP [XSt  1  ) .  ID<  1  ) 

REALM  RPERt  1  ) 

NEW  1 D  -  1 D <  1  > 

00  10  J-2.NUM 
OLD  1 0  -  IU< J  ) 

RPERt  NEW  1 0  )  -  RPER<  NEWID  )  ♦  RPERt  OLD  I  0  > 

RPERtOLDID  >  -  0. 

NP I XS ( NEW  I D )  -  NP I X  S  t NEWID  I  ♦  NP  I XS( OLD  ID > 

NPIXStOLUlD  >  -  0 

IF  t ISLSt NEWID >.GT. ISLStOLDID > >  ISLStNEWID)  -  ISLS(OLDID) 
ISLStOLDIOl  -  32000 

IF  <  ISSSt NEWID  )  .GT. ISSStOLDID > >  ISSStNEWID)  -  ISSStOLDIDl 
ISSStOLDID)  -  32000 

IF  <EL(NEWIDI.LT.EL(OLDID>)  E L < NEWID )  -  ELtOLDID) 

ELtOLDID  )  -  0 

IF  < ESt NEWIO  I  .LT .EStOLDID  )  )  ES(NEWID)  -  ES(OLDID) 

10  EStOLDID)  ■  0 

RETURN 
END 


FORTRAN  IV-P 
PCNCAT.FTN 


V02-5  1 

/TR : BLOCKS/VR 


19:92 'Si 


04-APR-85 


SUBROUTINE  PCNCAT< NCEN . PCNTR . BEG  1 N . END . NEW J D , RPE R  ) 

THIS  SUBROUTINE  CONCATINATES  ALL  ENDING  PERIMETER  VALUES  FOR  EACH 
I  .D.  ON  AN  ENDING  L INE . 

INTEGER  PCNTRI  1  ),|EGIN(  1  ),END<  1  I 
RFAL*4  RPER<  I  > 

DO  i 9ft  J-l .NCEN 

IF  <NCWID.EO.PCNTR(J M  GO  TO  101 
GO  TO  100 

H  RI’EIHNEWID)  ■  RPER<  NEWID  )  ♦  EHD<J>  -  BEGIN!  J  )  ♦  2. 

W  CUHTINUE 
RETURN 
END 


FORTRAN  IV-P 

FINO  FTN 


V02-51 

/TR : BLOCKS/WR 


19:03)00 


0 4-APR-05 


0001 


000 ?. 
00  03 
000  \ 
0005 
0000 
000  7 
0000 
000') 
00  10 
00  1  1 
00 1 ;! 

001  J 
nn  i  a 


SUBROUTINE  F I ND< NO  I N . NEWIO . NP I XS . OLD  I D  ) 
C 

C  USED  TO  FIMO  A  NEW  BIN  POSITION 
C  NBIN  -  LENGTH  OF  BIN  BUFFERS 
C  NEWIO  -  HEW  BIN  VALUE  RETURNED 

C  NP  I  XS  -  SUM  OF  PIXELS  BUFFER 

C  OLDIO  -  LAST  FOUND  I .0. 

INTEGER  OLDIO 
INTEGERM  HP IXS<  1  I 
DO  100  0-OLDID.NBIN 
IF  <  NP  I XS<  J  ) . EG .0  )  GO  TO  200 
100  CONTINUE 

DO  150  J-l. OLDIO 
IF  <NPIX3< J I.EQ.0)  GO  TO  200 
150  CONTINUE 

TYPE  *,  1  All  bins  f Iliad' 

STOP 

200  HEWIO-J 

OMJID-J 
RETURN 


FORTRAN  IV-P 
CORSRJ.FTN 
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/TR : BLOCKS/VR 


19:03:28 


04-APR-85 


1 


0001 


0002 
0003 
0001 
0005 
000C, 
000  7 

0000 
01)00 
00  10 
00  I  I 
00  1  2 
00  1  3 
00  1 1 
00  l'i 
00  10 
00  1  7 
0010 
00  19 
0020 
002  I 
0022 
0023 
0024 
0025 
00  20 
002  7 
002  H 
002  9 
0030 
003  I 
00J2 
0033 
0034 
0035 

0030 
003/ 
003U 
0039 
00  10 
001  1 
00  12 
004  1 
00  1 4 
00  15 
004  0 
00  4  7 
00  4  9 


P 


SUBROUTINE  CORSRT <KEY,PTR,LEN) 

C 

C  IN  CORE  SORT  ROUTINE.  SORTS  THE  VECTORS  HE V  AND  PTR  INTO  ASCENDING  ORDER  OF 
C  KEY.  THE  VALUES  IN  KCV  ARE  TREATED  AS  LOGICAL  QUANTITIES  HENCE  SION  BITS 
C  MUST  BE  TREATED  ACCORDINGLY.  ALPHABETIC  RECORDS  MAY  BE  SORTED  SINCE 
C  CHARACTER  CGUCS  ARE  IN  LOGICAL  OHULR.  THE  PTR  ARRAY  CAH  THEM  BE  USED  TO  MOV 
C  RECORDS  IN  A  OISK  F  ILL . 

C 

IMPLICIT  INTEGER! A-Z  > 

INTEGER *4  IMSK  .  JMSK . SM< 32  ) , PTR <  I  l.TEMP 
LOGICAL *4  I MSL  .  JHSl.  .KEY!  1  I  .L TEMP 
DIMENSION  BOR Y( 3 . 32  )  , COD ( 32  I 
EQUIVALENCE  !  IMSK.  IM'.iL  ).(  JMSK.JMSL  ) 

DATA  SM(  1  )  .  3M<  2  )  . SM<  1 7 1/020000000000,010000000000.0100000/ 

C 

00  1  1-3.16 

1  3M<  I  )  -  SM<  I  -  1  )/2 
DO  2  1-10.32 

2  SM (  I  )  -  SM<  1-1  1/2 
LEV  -  I 

DDK Y! 2. LEV)  -  1 
BDU Y <  J , L L V  I  -  LEM 
CUD! LEV  I  -  2 

72  CO  -  CIIOT  LEV  ) 

PL  -  BDR Y<  CD . LEV  I 
PH  -  BURY! CB+ 1 .LEV) 

IF  <PL.GE.FU>  GO  TO  75 
IMSK  -  3M<  LE V  ) 

01  JHSL  -  KEYf PL  > .AND . IMSL 

IF  <  JMSK.  HE.  0)  GO  T(J  US 
Pi.  -  PL  •»  I 

IF  (PL .CU.PU)  GO  TO  73 
GO  TO  81 

05  OHSL  -  KEY! PU). AMO. IMSL 

IF  (OMSK. HE. 0)  GO  TO  86 
LTEMP  -  KEY! PL  I 
TEMP  -  PTR! PE) 

KEY! PL  I  -  KEY! PU  ) 

P  TR( PE  )  -  P  TR  < FU  ) 

KEY(PU)  -  L TEMP 
PTR(FU)  -  TEMP 
GO  TO  81 

86  PU  •  PU  -  1 

IF  <PU.NE.Pl  )  CO  TO  05 

73  JH3E  -  KEY<PL  I. AND. IMSL 

IF  <  JMSK  .NE  .0  I  PI.  ■  PL  -  1 
IF  ! LEV  ,GE  .  32  )  GO  TO  75 
CB  -  CUD! LEV  I 
LEV  •  LEV  ♦  I 
CUD! LEV)  -  1 

BURY :  1  , LCV  )  •  BORYfCO. LEV-1 ) 

BURY! 2.LLV  )  -  PL 

IIDRV!  J.I.L'V)  -  BDRY(C0*1  .LEV-1  ) 

GO  TO  72 

75  IE  (CI)U<  LEV  >  .EQ.2  )  GO  TO  76 

CUD! L L V )  -  2 


-A  13- 


PORTRAIT  IV-P 
CORSR7.F7N 


vat-',  i 

/TR : ULOCKS/WR 


19:03 >28 


04-APR-85 


P 


ODRY<  2 . LEV )  -  BDRY( 2 .LEV  ) 

G U  70  72 

LEV  -  LEV  -  1 

IT  <  LEV .GT .0 )  GO  TO  75 

nr.  turn 

END 


FORTRAN" IV-P  V02-51  19i*l;54  04-APR-B5  P  1 

ENOFL  l.FTN,  /TR : BLOCKS/VR 

0001  SUBROUTINE  ENOFL  HRPER.NPIXS, 1SLS, ISSS.EL.ES.NL.LINE, 

1  OLCNTR  .  MECNTU ,  NOl.C  ,  HEWC  .  MUM  ,  OBEG  1  N  ,  NBEG  I M  . 

2  OEND.NEND.AUUF .PUUF .L  1M  IT  ) 

C 

C  SEARCH  FOR  FINISHED  OBJECTS 
C 

0002  1NTEGERM  NPIXS(l) 

000:1  INTEGER  OlCM7R<  1 ) . MECNTRI  1 ),OBEGIN( 1 1 ,NBEGIN(  1 ).OEND< 1 ) . NEND( 1 ) 

0004  INTEGER  I  SL S (  1  > , I 3SS<  1  ) .  CL (  1  )  ,  E S(  1  ) 

000'i  REALM  RPER<  1  ),ABUF<  1  >,PUUF(  1  ) 

000b  IHTEGER*2  L I H I T <  4  .  1  » 

C 

C  SEARCH  OLCNTR  FOR  I.U.'J  MOT  IN  NECNTR 
C 


0007 

IF  (LINE  .EG.  HI.  )  GO  TO  52 

00011 

IF  (NOLC.EO.0)  RETURN 

0009 

IF  (MEWC.EO.0)  GO  TO  51 

00  10 

l)U  10  J-l. NOl.C 

001  1 

NEW  IN  -  OLCNTR(J) 

0012 

UO  20  K-l.HEWC 

00  l  3 

IF  (NrWID.EU.HECNTR(im  GO  TO  10 

0014 

r 

20 

CONTINUE 

C 

p 

DROPLET  HOT  CONTINUED 

00  lr. 

CALL  PCNCAT < NULC, OLCNTR, OBEG1N.OE NO, NEWI D , RPFR ) 

00  1  b 

CALL  F  I  HAL  <  HEW  1 1) ,  RHER  ,  NP  I XS  .  I  SL 3  .  I SSS  .  EL  .  E 3  .  HUM  ,  ADUF  . 

1  PBUF. LIMIT) 

0017 

10 

CONTINUE 

0018 

r 

RETURN 

V# 

c 

p 

FIN] 

3H  OFF  EVERYTHING 

0019 

51 

DO  50  J-l.NOLC 

0020 

NEW  ID  -  OLCNTR! J ) 

002  1 

CALL  PCNCAT (NOLC, OLCNTR, OBEG I N.OEND. NEWI D. RPF R  ) 

0022 

CALL  F  INAL  <  NEWID  ,  III’ER , NP  I XS  ,  I3LS.  I S3S  .  EL  .  E3  ,  HUM  .  ABUF  . 

1  PBUF, LIMIT) 

0023 

50 

CONTINUE 

002  4 

RETURN 

0025 

52 

IK  (HOLC.EQ.0)  GO  TO  50 

002  ft 

IF  (HEWC. EG. 01  GO  TO  51 

0027 

DO  110  J-l.NOLC 

002U 

NEWID  -  OLCNTR ( J ) 

00  2*1 

DO  90  K- I , NE WO 

0010 

IF  (NEWID. EU. RECNTR(K) >  CO  TO  100 

003  1 

90 

CONTINUE 

0012 

CALL  PCNCAT! NOLC, OLCNTR. OBEG I N.OEND. NEW  ID, RPFR ) 

00  33 

CALL  F I NAL ( NEWI 0 . RPER , NP I XS , I3LS, I SSS . E L . E 3 . HUM .ABUF , 

1  PBUF. LIMIT) 

01)11 

GO  TO  00 

00  IS 

100 

CALL  PCNCAT! NEWC .NECNTR , NBEG IN. NEND. HEW  ID. RPER ) 

003  b 

CALL  F 1 NAL ( NFWI U. RPER. NPJXS, I SLS . I 333 , E L . L 3 . HUM . ABUF . 

1  PUUF. LIMIT) 

0017 

00 

COMTIHUC 

00  30 

RETURN 

rO*T*AN“IV-P 

ENOFLt.FTN 


V02-51 

/TR  : ULOCKS/WR 


19:01:54 


04-APR-85 


0039 
niuts 
00  4  i 
0042 
004  3 

004  4 
00  4'i 
0045 


P 


60  IF  <  NEVC . EO . 0  )  RETURN 

1)0  70  j-i,mi:wc 
Nil VI U  -  MECNTRI J  ) 

CALL  PCNCAT  <  NEVC  ,  MECNTR  .  NflEG  I N  ,  NE NO  ,  NEW  t  D  ,  RPF R  > 

CALL  F  INAL  I  MEV1D  .  Rl’CR  .  Nl’  I XS  ,  1SLS.  1SSS  .  EL  .  E S  .  HUM .  ABUF  , 
1  PBUF. LIMIT) 

70  CONTINUE 

III  TURN 
EMU 


19:29*49 


25-MAY-82 


PACE  1 


FORTRAN  IV-PUJS  V02-51 
FINAL. FTN  /TR:BLOCKS/VR 


0001 


nr*  at 

iii  '»:i 
min  < 
ni'ii'i 
jl'iilli 
iinn  ' 
mill 
ii'Hi'i 

mn  n 
mi  1 1 
DU  1 2 


nr m3 
mi  1 1 
lllll’i 
HU  1  li 
mi  i ; 

IllJ  1  u 

ll'll'l 
1(020 
nu2  i 

0022 

IMI2T 

Hi*24 

i)"2'i 

I1D2I) 

i)i).’.  / 

uii2n 

1(1)2’) 


SUOROUT 1NE  F INALINFVIO.RPER.NP IXS. ISIS, I SSS . EL . E S . MUM  .  AQUF  . 

I  PUUF. LIMIT) 

C 

C  TO  TERMINATE  AN  OBJECT 
C 

INTEGER *4  NP  IXS(  I  ) 

I MTEGE R*2  LIMITS.  1  > 

INTEGEK  ISLSI  1 >, I3SSI I  )  ,EL( 1 >.E3(  1  ) 

RCAI.  *  4  RPER1  1  ).AUUF<  1  )  , PUUF 1  1  ) 

U  VTL  OR l GA 1 5 1 2  > 

COMMON  /MI HARE/  AMINAR 

II  < HP  1X31 HEVIU ) .EO.0)  RETURN 
HUM  -  NUM  ♦  1 

IF  <NUH.LT.2JM  )  GO  TO  5 

TYPE  MuinUor  of  spots  gt  buffer  space  available' 

S  f OH 
C 

C  ACCUMULATE  DROPLET  STATISTICS 

C 

5  ARUF ( HUH  )  -  NPIXSINCVID) 

HI1UF (  HUM  )  -  RPER1  NEW II) ) 

L  I M  I T  (  1  .NUM)  -  131.31  NEVID  I 
LJMrr<2,NUM>  -  IS3S1  MI'VID  ) 

L  IM  IT! 3  ,  HUM  )  -  El.  1  NEW  I U  )  -  tSLS<NEVID)  ♦  1 
I.  IM  IT!  4  .  HUM  )  -  ESI  NEW  ID  >  -  1SSS1HLVID)  ♦  1 
C  CLEAR  OUT  H I  M3 

IE  < ARUF < MUM >.GT. AMINAR)  GO  TO  10 

CALL  IIORKARl  I  .  L  IM  I  T  <  1  .  MUM  ) .  L  IM  I T  (  2 .  NUM  > .  L  IM  IT<  3  .  NUH  ) .  L  IH I  T<  4  .  NUM  M 
HUM-HUM- I 
C 

10  HP  1 X3 1  MEW  I D )  -  0 

IIPEKI  HEWiO  )  -  0. 

I  31.31  HEVIU  )  -  32000 
1 3331 NEWID  )  •  32000 

El.  1  MKV1  ,j  )  -  0 

ESI NEWID)  •  0 

RE  ruiiH 

EMU 


-A17- 


W  v  v  v '  j'.v,  /  i ,-.v  -V-V-  /  j .p.v.  •fw-f. 1 -r.  ■>'. 


FONTMNMV^P 

MINST. FTN 


V02-51 

/TR:BLOCKS/VR 
PROGRAM  RA  INST 


10:50:30 


23-FEB-04 


Road  droplots'  data  file  and  compute  general  statistic 

Programor:  Miguel  A.  Ilornan 

REAL  * 4  nniJF<  1200)  .  V0I1UF  <  (200>.CSFFR?  1200) 

It  I.AE  *0  '•HU') .  SH02  .  D TO  f .  E3T0 T  .  5HU4  .  C0H0U 

imteci:r*2  muhc 

LOG  I CAL  *  1  AN  SA 

REAL -4  VOI.II  1 3?  0:  79  )  .FREI1  IS?  0:  79  ) 

REAL*  4  AREA  .  PER  IH  .  1)3 . 02 . 1) . D4 

l.o*;  I  r  AL  *  1  HO.'OIM  12  >  .MO2IM?(S0  > .NOZID?  00) 

ISOM  (VALENCE  <  M020IK  1  ).M021U<  1  )  ) 

EDI)  (VALENCES  N021M?  l  ),M0210<  13)  ) 

IN  TECER-2  OPE  I'Ll) 

I.OOICAI.M  F  1 1. E HA <  22  > 

INTEGERS  1MUFCK?  H )  . SI  . 53 , ML  . NS 
E*)M  (VALENCE'.  IIIOFCIl?  I  )  ,SI.  ) 

coo  i  vale  nee  ?  iiiuecu?  2  > .  v; ) 

CUD  I  VALE  ME  El  IIIUCCK?  I  )  .Ml.  ) 

E  OH  (VALENCE?  IIIDFEK?  1  >  ,M3  > 

COO  I  VALENCE?  I  HUE  Kit?  ’>  I  .AREA  ) 

EOU  (VALENCE?  IIIHFEIt?  /I.PERIH) 

DA  ?A  VOI.II  13/00*0.  /  .  FKE1U3/I) 0*0.  / 

DATA  N020D  /  1  N ’ , ’ 0 ' . ' 2 ’ . ' 2 ’ . ’ L ‘ ’ E ’ . '  I  D  '/ 

Open  data  file 
FORMAT (0.22.M  ) 

TYPE  •«’  Enter  filename' 

ACCEPT  45 1 .  LONGU . <  F I  LENA? I  > . I- 1 ,LOHGO ) 

OPEN? UNIT-3, NAME «F  ILENA. TYPE- 'OLD'  . 

0)1  SP"  KEEP  '  . ACCESS- 'DIRECT ' , RECORDS 1ZE-4 . FORM- ' UNFORMATTED '  ) 
READ?  3 ' 1  )  II1UI  ER 
NUHE-IIUJFER?  1  ) 

ICON? -0 
•»M'M=0. 

SMI)  1-0. 

SMD2-0. 

I)  TOT-0. 

E STOT-0. 

PI-:).  14159 

01-4 ./<  3. -SORT?  F I ) > 

Main  loop 
DO  100  K-2.HUME 
READ?  3 'K)  IllUFF.R 
II  ?  IIIUFEN?  1  I.EU.0)  GO  TO  200 
ICONT «  ICONT ■»  I 
D2-? 4 . -AREA  )/P  I 
I)  *SOR  f  <  1)2  ) 

D  I  0)2*1) 

I)  I -112-1)2 

E 'SEE. It?  1C0MT  )-<PCRIM*PERIH)/<AREA*4.-PI  1 
Dll  HE  ?  (COMT l-D 


FdkTRAN 

"1V-P 

V03-ril  10:50130  23-FEB-04 

RAlNST. 

f  TN 

/ f K i ULOCKS/WR 

mur, 

VOIHJF  <  ICONT  I »Q I  * ( ARC A**  1 .5  1 

')•)  1  / 

shdi-smim+im 

l*i*  in 

;iii)J-:;hi)J-mu 

o<t  i '* 

SHD2-3MI)2  +  I)2 

iliVid 

i)  ro  r-uTOT-M) 

lilts  1 

es tot-es  ro  t+esfcrt  icont  i 

I»«1S? 

1  1-10*1) 

(Ill'll 

I  (-•MINOT  7'i.ll  > 

ll'I'i  I 

i  uiui i •;<  1 1  )=*ruciiis<  !  I  >-*i 

I'O’i'i 

VOI.IIIST  1  I  I-VOLIIIST  I  I  >-*03 

tlil'ili 

1  00 

con  r  i  not: 

ll'I'i/ 

CLOSE  (UNIT-3) 

1 1 1  * '  i  •'  1 

NO  TO  300 

llfl'i'l 

200 

CLOSE? DM  IT-1 > 

1)060 

TYPE  *.'  Kilo  records  counter  error ' 

1)1)0  1 

f 

CALL  EXIT 

l, 

i: 

Uegln  general  statistic 

noc? 

i, 

300 

TYPE  *,’  Enter  nozzle  coda' 

Ill) li  » 

ACCEPT  750,  ITEX.IHOZIMT  1  >.  1-1.  HEX) 

ll'I'i  t 

750 

FORMAT (  '.lillAl  > 

It'lli’i 

TYPE  *,'  Enter  watar  pressure  <pal)' 

ll'lii'i 

ACCEPT  * , WATP 

ll'I'i/ 
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